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^ ! 

Q ■ I give a brief overview of our present understanding of the high temperature phase of 

^ . QCD, trying to clarify some of the theoretical issues involved in the current discussions 

^ ! that emphasize the strongly coupled character of the quark-gluon plasma produced at 

> • RHIC. 

^ : 
o . 

I ■ 1. From the "ideal gas" to the "perfect liquid" 

^ ■ QCD asymptotic freedom leads us to expect that matter at very high temperature 

Q . and/or density should become simple [TP] : an ideal gas of quarks and gluons whose free 
"rH ! motion is only weakly perturbed by their interactions. Lattice calculations indeed show 
! that, at very high temperature, the thermodynamical functions indeed go, albeit slowly, 
Cj I towards those of free massless particles (see for instance [3]). We also start to have a 
^ ■ semi-analytic understanding of this high temperature state in terms of quasiparticles, 
^ . as we shall discuss later. In agreement with general theoretical considerations on the 
• ^ ! running of the effective coupling with the temperature, lattice calculations also show that 
^ I as the temperature decreases, the coupling increases and the physics becomes less and 
' less perturbative. Eventually at some critical temperature, a phase transition, or a rapid 
cross-over, takes place towards a regime dominated by hadronic degrees of freedom. 

This standard, and theoretically well established, picture has been recently challenged, 
with most of the present discussions emphasizing the strongly coupled character of the 
quark-gluon plasma. Three elements have conspired to this shift in paradigm. First, the 
RHIC data do not provide any evidence for ideal gas behavior, but are better interpreted 
by assuming that the produced matter behaves as a liquid. Second, perturbation theory is 
notoriously unable to describe the quark-gluon plasma unless the temperature is extremely 
high. Third, new techniques have emerged that allow strong coupling calculations to be 
done in non abelian gauge theories. 

Most proposed "signatures" that have been suggested during the preparation of the 
heavy ion program, are based on the simple picture of the quark-gluon plasma as an ideal 
gas of weakly interacting constituents. Part of the reason for this is simple: this is a 
regime where one can do first principle (and elementary) calculations. However, RHIC is 
forcing us to look into a region where theory is hard, i.e., study the quark-gluon plasma 
in a regime where the coupling is not very small. The temperatures reached at RHIC are 
presumably not high enough to reach the asymptotically free regime. We learned from 
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RHIC (see [illSHGlTf] . and also [8]|9] ) that the matter produced in ultrarelativistic heavy 
ion reactions has a large energy density, high enough for the early stages of the collisions 
to be dominated by partonic degrees of freedom. Suppression of jets reveals also strong 
densities and strong interactions. The observed elliptic flow indicates that the matter 
behaves collectively, as a liquid with low viscosity [lOj. Although it is difficult to place 
quantitative bounds on the values of the viscosity from the analysis of the data, it seems 
indeed that low values of the viscosity are needed PJJ, values that seem to be incompatible 
with perturbative QCD calculations. 

Much effort has been put into calculating the successive orders of the perturbative 
expansion for the pressure p^l3lll4p5] and the series is known now up to order In g [15] . 
Such calculations indicate that the asymptotic quark-gluon plasma indeed emerges, as 
expected, at very high temperature. However, they also show that perturbation theory 
makes sense only for very small couplings, corresponding to extremely large values of T . 
For not too small couplings, the successive terms in the expansion oscillate wildly and 
the dependence of the results on the renormalization scale keeps increasing order after 
order [16]. This situation is to be contrasted with what happens at zero temperature, 
where perturbative calculations achieve a reasonable accuracy already at the GeV scale. 
At finite temperature, thermal fiuctuations alter the infrared behavior in a profound way. 
Nevertheless, as I shall argue, this does not imply that weak coupling calculations are 
useless. In fact a lot can be learned from them, in particular how to capture the right 
physics that can allow for smooth extrapolations to the strong coupling regime. 

Finally, new techniques for doing calculations in strong coupling have emerged recently, 
giving hope that we may soon have at our disposal new tools to handle strongly coupled 
gauge theories. These techniques are based on the recognition that some supersymmetric 
Yang-Mills theories are dual to gravitation theories. The duality involves an interchange of 
the regimes of weak and strong coupling: weak coupling in the gravity theory corresponds 
to strong coupling in the gauge theory. This duality offers the possibility to study strongly 
coupled gauge theories by performing perturbative calculations in their gravity duals. This 
possibility has been exploited in a number of recent publications. One prediction of such 
calculation is that the entropy S behaves, in strong coupling as [17] 

-5 3 45,, , 1 

where A = 2g'^Nc and 5*0 is the entropy of the non interacting system . Thus, in the limit 
of strong coupling, A — > cxo, the entropy is bounded from below by the value S/ Sq = 3/4. 
The fact that this value is close to that obtained in lattice calculations at temperature 
above Tc has contributed to sustain the speculations that the quark-gluon plasma above 
Tc is in a strongly coupled regime. Note however that at T ~ ST^, the value of the 
entropy density is half-way between its weak coupling value and its strong coupling value 
(see Fig. [1] below), so there is no compelling reason to favor at that point an approach 
based on a strong coupling expansion. Furthermore, the minimal value of 3/4 is obviously 
not compatible with lattice data near where the entropy vanishes and the coupling is 
the strongest. In fact, supersymmetric Yang-Mills theories have symmetries that make 
them rather different from QCD: in particular the coupling constant does not run, and 
there is no phase transition. The running of the coupling constant is however an essential 
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property of QCD; in particular, above the transition region it is accompanied by a breaking 
of conformal symmetry that is very well observed in lattice calculations P^19f2U] . 

In fact the physics above T^, in the temperature range ^ T ^ 3Tc, remains poorly 
understood, with many unanswered questions: What are the degrees of freedom relevant 
for an effective theory? Are remnants of confinement important, such as those captured 
by effective theories for the Polyakov loop (see for instance [2T||22] or, on a more phe- 
nomenological side, [23])? What is the role, if any, of bound states [21]? What is their 
fate as the temperature grows, do they survive at large temperature, as recent lattice data 
suggest [25]? What are the charge carriers [261127] ? Etc. 

In the rest of this talk, I shall not discuss this transition region, but focus on the high 
temperature phase. Also I shall only consider the regime of vanishing chemical potential 
(for a recent study of QCD thermodynamics at finite chemical potential, see [28]). 

2. Weakly or strongly coupled quark-gluon plasmas 

In order to get orientation into the effects of the interactions a general strategy is to 
compare the kinetic energy of the particles to their interaction energy. When the kinetic 
energy dominates, the non-interacting system constitutes a good starting point for the 
theoretical description, and one may think of using perturbation theory. Note that this 
kind of argument may be obscured by the emergence of new degrees of freedom in the 
interacting system: in this case one could be led to conclude that it is strongly coupled, 
while the main effect of the interactions is to generate these new degrees of freedom 
which may have weak residual interactions. The natural starting point of the theoretical 
description in such cases would be the non-interacting system built with these new degrees 
of freedom. 

Consider an elementary plasma, made of electrons, positrons and photons. When the 
charged particles behave individually as classical particles, the mean kinetic energy of an 
electron is simply proportional to the temperature T, while the average potential energy 
per particle is of the order e^n^/^, where n is the number density, so that n^^^ is the inverse 
of the interparticle distance. The condition that the plasma be ideal is then T 3> e^n^^^. 
The dimensionless parameter = e^m}/^ /T is essentially the plasma parameter F used 
in plasma physics [2^. In terms of the condition that the plasma is ideal is simply 
5^ -C 1. The Debye screening length is A ~ \Jt /ne^ ~ n~^/^ / g. Thus when g is small, the 
screening length is large compared to the interparticle distance, which is a criterion for 
collective behaviour. Also when g is small individual collisions can be ignored. 

In ultrarelativistic plasmas, the temperature and the density are no longer independent 
control parameters since ~ T'^. Then the plasma is characterized by a single dimension- 
full parameter, T. This is the situation in QCD at high temperature, where the parameter 
g is the gauge coupling (at a scale ~ T). Note that, just above Tc, the couping constant is 
not small, but not huge either, (7 ~ 2 (see e.g. ^30]). However, to decide whether the quark- 
gluon plasma is strongly or weakly coupled it is essential to recognize that the effects of the 
interactions depend not only on the strength of the coupling but also on the magnitude 
of the thermal fluctuations, which depends on their wavelengths. At weak coupling a 
hierarchy of scales of thermal fluctuations emerges, each scale being associated with well 
identified physics [31] • Predicting the effect of the interactions amounts to comparing 
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the kinetic energy ~ {{dA)"^) ~ k'^{A'^) with interaction energy g'^{A'^) ~ g'^{A'^)'^, or 
equivalently k"^ with g^{A'^). 

The plasma particles have typical energy and momentum of the order of the temperature 
k ^ T. The thermal fluctuations at this scale are (^4^)7- ~ T^. These hard fluctuations 
constitute the dominant contribution to the energy density, but they produce only a 
small perturbation on the motion of a plasma particles if g is small. Consider now a soft 
excitation at the momentum scale k ~ gT. The fluctuations at this scale are {A'^)gT ~ 
gT^. The self-interactions of soft modes remain a small correction {g'^{A^)gT ~ g^T"^ <C 
/c^ ~ g'^T'^), but the interaction of the soft mode with the hard fluctuations cannot be 
ignored {g'^{A'^)T ~ g'^T'^ ~ k"^)'- the hard particles renormalize non-perturbatively the 
propagation of the soft modes; this is the physical origin of the hard thermal loops. And 
indeed, the scale gT is that at which collective phenomena develop. The emergence of 
the Debye screening mass m^) ~ gT is the simplest example of such phenomena. Moving 
down to a lower momentum scale, one meets the contribution of the unscreened magnetic 
fluctuations which play a dominant role for k ~ g^T. At that scale, {A^)g2x ~ g'^T'^i 
so that g'^{A'^)g2rp ~ g'^T'^ ~ k"^: the fluctuations are no longer perturbative. This is the 
origin of the breakdown of perturbation theory that occurs at order g^ in the pressure 

This particular pattern of scales, and the way they are coupled, is what complicates 
the theoretical description of the quark-gluon plasma. Even when the coupling is weak, 
non-perturbative effects arise: for instance collective modes develop, whose treatment 
requires resummations or the use of effective theories. The real issue here is not so much 
the strength of the coupling but the fact that there are many quasi degenerate active 
degrees of freedom. As we shall see in the rest of the talk, weak coupling techniques allow 
us to identify and perform the appropriate reorganizations of the perturbative expansion 
that are needed in order to take into account the role of the new degrees of freedom, in 
particular the collective modes. We shall see that once this is done, the extrapolation to 
strong coupling is much smoother than what perturbation theory could lead us to expect. 

3. Weak coupling techniques for strong coupling regimes, resummations 

Since the major problems in thermal QCD come from the infrared sector, and in partic- 
ular from modes carrying zero Matsubara frequencies, it is natural to isolate those modes 
and construct for them an effective theory through a procedure known as dimensional 
reduction. The resulting Lagrangian is of the form [32|33|34j 

£e = ^ tr 4 + tr [A, Ao]^ + m| tr Al + ^Ai,( tr Alf + ... (2) 

where the parameters gE {Dj = dj + igEAj), uie and A^; are determined perturbatively as 
a function of the gauge coupling g by matching. In lowest order g'^ = g^T and m|; ~ g^T^, 
Xe ^ g^T. 

Calculations based on this scheme have been pushed to high order [35] , but they depend 
on an as yet undetermined 4- loop matching coefficient. By adding a parameter to account 
for the uncalculated g^ contribution, one can match the four-dimensional lattice results. 
The required value of the coefficient is not very large, suggesting that the contribution of 
the magnetic sector to the pressure is numerically not important at high temperature. 
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Figure 1. Comparison of the lattice data for the entropy of pure-glue SU(3) gauge theory 
of Ref. [TSj (gray band) with the range of Shtl (solid lines) and S^la (dash-dotted lines) 
(see text). From Ref. [32|. 




An interesting feature of the perturbative calculation when organized through the di- 
mensionally reduced effective theory ([2]) is that the large scale dependences of strict per- 
turbation theory can be reduced when the effective parameters are not subsequently ex- 
panded out [T5|36] . with noticeable improvements when going from two-loop to three-loop 
in the effective theory. 

Other ways to reorganize the perturbative expansion have been tried. One proposal, 
called "screened perturbation theory" [371138] . has been generalized to QCD by Andersen, 
Braaten, and Strickland A more systematic approach to screened perturbation 

theory is based on an expression for the entropy density that can be obtained from a 
^-derivable two-loop approximation [^TP2] . This approach has a clear physical content: 
the dominant effect of the interactions is to turn the original degrees of freedom, quarks 
and gluons, into massive quasiparticles, with weak residual interactions. This theoretical 
approach gives support to the more phenomenological one proposed in Ref. [13] 

In Refs. pTIH2] the lattice results for the thermodynamic potential for T > ST^ were 
quite well reproduced as can be seen in Fig. [H In this figure, numerical results are given 
for the entropy obtained in two successive approximations: in the first, denoted HTL, 
the hard thermal loop self-energies are used as an input of the calculation, in the second, 
denoted as NLA, corrections to the HTL self-energies are taken into account. The full 
lines show the range of results for Shtl when the renormalization scale p, is varied from 
ttT to 4:7iT; the dash-dotted lines mark the corresponding results for Snla- The dark-gray 
band are lattice data from Ref. [18] (the more recent results from Ref. [H] are consistent 
with the former within error bars). Evidently, there is very good agreement for T ^ 3Tc. 
The result from Ref. [SUPU] is indicated by the dooted lines at the top of Fig. [TJ 

The 2PI formalism used in the calculation of the entropy that we have just described 
has been tested 05] in a model with a large number of quark flavors, which can be 
solved exactly. In this case, the 2PI formalism that we used is exact. The comparison 
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provides then indications on the quahty of the further approximations involved in solving 
the 2PI equations. The results are quite encouraging, and a remarkable agreement can 
be achieved for quite large values of the coupling constant. Note however, that these 
calculations indicate that much accuracy can be gained by taking fully into account the 
momentum dependence of the corrections to the asymptotic thermal masses that enter 
the calculation of the NLA (in the results reported in Fig. [1], only momentum averages 
have been used [11112] )■ Implementing such corrections is technically demanding and has 
not been completed yet. 

4. Insights from the functional renormahzation group 

In the last part of this talk, I would like to examine the issues addressed above from 
the perspective of the functional renormahzation group (fRG) [I5|47|48|49II5U] (sometimes 
called exact or non-perturbative, depending on which aspects of the formalism one wants 
to emphasize). 

The fRG has been applied, in various incarnations, to a variety of problems (for reviews 
see [?T|52] ). It allows for non-perturbative approximation schemes which can, in some 
instances, be quite accurate. Functional RG techniques have been applied to problems 
at finite temperature in the past P8|53|54|55|56j . and the general behaviors that we shall 
discuss here have been known already for some time. However, the main focus of previous 
studies has been the description of the phase transitions, rather than the specific problem 
that we address here. 
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Figure 2. The flow of the coupling constant as the function of the ratio k/A, where n is the 
scale of the infrared cut-off, while A is the microscopic scale. The left figure corresponds 
to weak coupling {g{2TxT) = 0.0994), while the right figure corresponds to strong coupling 
{g{2nT) = 1.1). 

There is some analogy between the effective field theory approach and the functional 
renormahzation group: in effective field theory one integrates out degrees of freedom 
above some cut-off; in the renormahzation group this integration is done smoothly. The 
general strategy of the fRG is to control the long wavelength fluctuations with an infrared 
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regulator depending on a continuous parameter k. When k is of the order of the micro- 
scopic scale, fluctuations are essentially suppressed. On the other hand when k — they 
are all included. The presence of the regulator allows us to write flow equations that can 
be integrated all the way down to /t = in order to get the physical quantities. In a 
way, the renormalization group builds up a continuous tower of effective theories that lie 
infinitesimally close to each other and are labeled by the momentum cut-off scale k. These 
effective theories are related by a renormalization group flow equation. This picture is 
independent of the value of the coupling, and as we shall see, the renormalization group 
provides a smooth extrapolation from the regime of weak coupling, characterized by a 
clean separation of scales, towards the strong coupling regime where all scales get mixed. 

I shall only present a few results obtained recently together with A. Ipp, R. Mendez- 
Galain, and N. Wschebor [57]. These results concern a scalar (j)'^ theory with 0{N) 
symmetry which exhibits similar bad convergence properties of its thermal perturbative 
expansion as QCD. In our study, we have used the so-called local potential approximation 
(LPA) [52j. We have shown that this approximation is compatible with perturbation 
theory up to, and including, order g^, and that is provides a smooth extrapolation to 
strong coupling similar of that of a simple 2PI approximation. (Note that in scalar 
theories, the presence of the Landau pole puts a limit on the maximum value of the 
coupling for wich calculations are meaningful. Strong coupling here means g ^ 1.) 
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Figure 3. Pressure as a function of the coupling g{27TT). The result of the RG-LPA 
calculation is compared with the perturbative estimates of order g"^ and g^ and also with 
a simple 2PI resummation. 

The left panel of Fig. [2] illustrates the flow of the effective coupling as the momentum 
cut-off runs down to zero. One observes that as long as k ^ 27rT, the running of the 
coupling is not modified by the thermal fluctuations. These fluctuations start to play a 
role when k ^ 27rT, turning the flow into a three dimensional one. This provides, in this 
context, a nice illustration of the phenomenon of dimensional reduction. When k becomes 
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of order gT, the thermal mass provides an infrared cut-off and the flow stops (the infrared 
regulator playing then no role anymore). As revealed by the right panel of Fig. [21 when 
the coupling increases, the pattern does not change much: the thermal fluctuations start 
to play a role when k ^ 2nT, turning the flow into a three dimensional one. The range 
of values of k where the flow is three dimensional is reduced because the thermal mass 
increases with the coupling. At the same time the amplitude of the three-dimensional 
flow is enhanced by the larger value of the coupling. These competing effects contribute 
to the stability of the results as one goes from weak to strong coupling. 

This is reflected in the plot of the pressure as a function of the coupling, Fig. [31 The 
pressure calculated with the fRG is compared with the order g"^ and g^ of perturbation 
theory, as well as with the result of a simple self-consistent 2PI approximation (this 2PI 
approximation, as well as the approximate solution of the flow equations become exact in 
the large limit). 
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